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Scenarios for particle production in the GeV and TeV regions are reviewed. The expected in- 
crease with cm. energy of the average number of clans for the soft component and the decrease 
for the semihard one indicate possible classical and quantum behaviour of gluons, respectively. 
Clan thermodynamics, discussed in the paper, appears as the natural framework to deal with such 
phenomena. 
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I. INTRODUCTION 

Phenomenological analysis of multiparticle production 
in hadron-hadron collisions in the GeV region |jl|, || re- 
vealed interesting substructures, i.e. soft and semihard 
events, without and with minijets respectively, each class 
of events being described by a negative binomial (Pascal) 
multiplicity distribution (NB MD) with different values 
of the parameters, the average charged multiplicity n and 
k = / [D^ — n). D is here the dispersion. The weighted 
composition of the two MD's leads to the observed final 
charged particle MD. Clan structure analysis in terms of 
the average number of clans, N, and the average number 
of particle per clan, fie, allows to interpret nicely the on- 
set of above mentioned substructures in the GeV region. 

By extrapolating n and k behaviour to the TeV region, 
three possible scenarios [|l|, |^ for the semihard component 
have been investigated. The first one assumes that KNO 
scaling is satisfied both in the soft and semihard com- 
ponent. This situation should be compared with what is 
assumed in the other two scenarios for the semihard com- 
ponent where KNO scaling is strongly violated or has a 
QCD inspired behaviour through the cm. energy depen- 
dence of the corresponding NB MD parameters. Since, 
at the present stage of QCD, calculations of MD's and 
correlations in the GeV and TeV regions cannot be per- 
formed in a sound way, we can only rely on QCD inspired 
extrapolations of the parameters. The last two scenarios 
for the semihard component (the soft component is taken 
to be the same in all three scenarios), although more re- 
alistic than the first one, lead to a decreasing average 
number of clans and to the corresponding increase of the 
average number of particles per clan as the cm. energy 
increases. Since clans are independently produced by as- 
sumption it would be important to understand the real 
meaning of their decrease for cm. energies in the TeV 
region, a fact which seems to widen the motivations at 
the origin of the first introduction of clan concept in high 
energy phenomenology ^, |^ . 

In addition, it should be pointed out that clan struc- 
ture analysis can be generalised to the huge class of dis- 
crete infinitely divisible distributions (IDD), to which 
NB(Pascal)MD belongs, and therefore any result ob- 
tained in the framework of clan structure analysis can 
be easily extended to the full class of IDD's. 

In the present paper after a short introduction on the 



most appealing statistical theories of multiparticle pro- 
duction a new interpretation of n charged particles multi- 
plicity distribution, p{n), for the class of IDD is proposed 
in terms of the canonical and grand-canonical partition 
functions. Then the connection of this new interpretation 
for the description of the soft and semihard substructures 
in hh collisions in terms of the NB (Pascal) MD is exam- 
ined. 



II. AN ALTERNATIVE APPROACH TO 
MULTIPARTICLE PRODUCTION 

One of the best known statistical approaches to multi- 
particle dynamics is Feynman's fluid analogy ^ , where 
the cross-section for the production of n particles plays 
the role of the partition function in the canonical ensem- 
ble, as it is an integral over phase-space of the square of a 
matrix clement which plays the role of the Gibbs distribu- 
tion, g-H/ksT ^ -j-j^jg approach the volume is identified 
with the extension of phase space and the fugacity z with 
the dummy variable u appearing in the definition of the 
generating function G: 



G(u)^^^.Xn). 



(1) 



This identification is unsatisfactory because one has to 
satisfy at the same time the definitions of the average 
number of particles, from the grand-canonical ensemble: 



in) 



dlnG 
dz 



and from the definition of generating function: 



dG 



rflnG 



u=l 



du 



(2) 



(3) 



The above formulae can be satisfied at the same time 
only in the limit of zero chemical potential. 

Another approach was proposed by Scalapino and 
Sugar [0 : they defined the probability amplitude to pro- 
duce a particle at rapidity y, denoted by n(j/), as a ran- 
dom field variable, then introduced a functional F[Il] 
which played a role analogous to the free energy for a 
system in thermal equilibrium. One can then obtain the 
n-particle inclusive distribution by averaging the prod- 
uct of the squares of the amplitudes, W^{yi) • ■ •n^(j/„). 
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with a weight given by e~^^^\ Lacking the knowledge 
necessary to calculate i^[n] from the underlying dynam- 
ics, the authors parametrised it (following Ginzburg and 
Landau) in retaining the first three terms in a series ex- 
pansion, then solved the model in a few particular cases. 
Remarkably, to leading order in the size of the allowed full 
rapidity range, they obtain a generating function which 
has the form of an IDD. 

More recent results obtained in the above mentioned 
frameworks, concerning KNO scaling and phase transi- 
tions, can be found in |, |, |lO[ ^ ^ |l| . 

Stimulated by these results we propose a new simpli- 
fied approach to the statistical theory of multiparticle 
production, heavily based on IDD properties and valid 
for any chemical potential. 

We denote with Qn{V, T) the partition function in the 
canonical ensemble for a system with a fixed number 
n of particles, volume V and temperature T, and with 
2(0, y,T) the grand-canonical partition function for a 
system with fugacity z, volume V and temperature T; 
the chemical potential ^ is defined by z = eicp{fj,/kBT), 
where fc^ is Boltzmann constant. 

We recall the relation between the partition functions: 



Q(z,l^,T) = ^z"g„(l/,T). 



(4) 



n=0 



Accordingly, in a statistical mechanics framework, the 
probability p{n) of finding n particles in the system is 
the following: 



p{n) 



Q{z,V,T) 



1 z" (9"Q 



(5) 



z=0 



Noticing that Qo{V,T) = 1, we find immediately that 
the grand-canonical partition function is the inverse of 
the void probability p{0), i.e., of the probability to find 
no particles in the system: 



p{0) = [Q{z,V,T)]- 



(6) 



This result is very general: the void probability is the 
inverse of the grand-canonical partition function and all 
properties of the system can be obtained from it 1 1 



Consider now the wide class of power series distribu- 
tions (PSD), usually defined as follows: 



p{n) 



7(6)' 



(7) 



with a, 

tion requires that 



and b free parameters, while proper normalisa- 

lirpw that 



7(&) = ^a„&". 



(8) 



n=0 



Notice that ao can always be chosen to be 1 (by redefining 
7(6) as 7(6)/ao). Then one has for the void probability: 



p(0) = 



lib)' 



(9) 



Comparing Eq.s (||) to (|^), our new approach is char- 
acterised by the following correspondence: 



(10) 



i(b)^p{oy 



A very interesting property of this novel identification 
is that a„ is the canonical partition function for a system 
with a fixed number of particles n, and in particular ai is 
the canonical partition function for a system with 1 parti- 
cle. This means that if we know the multiplicity distribu- 
tion of a thermodynamical system, and cast it into a PSD 
form, wc can not only deduce the grand-partition func- 
tion but also identify the fugacity of the system and the 
canonical partition function. As an intriguing example of 
this correspondence, motivated by the phenomenological 
analysis of multiparticle production in the GeV region, 
we will in the next section examine the NB(Pascal)MD. 



III. THE NEGATIVE BINOMIAL 
DISTRIBUTION 

Any discrete infinitely divisible distribution (IDD) can 
be written as a compound Poisson distribution (CPD) 
]l6t , i.e., the number of clans N can be defined in such a 
way that the void probability is: 

p{0) = exp(-7V). (11) 

Comparing with Eq. (||), we notice that for any IDD the 
average number of clans is the logarithm of the grand 
partition function: 



N ^\nQ. 



(12) 



All thermodynamical properties can then be obtained by 
differentiating the average number of clans. In particular, 
being for the grand canonical ensemble PV = ksTln Q, 
we obtain the following equation of state: 



PV = NksT; 



(13) 



it tells us that our system behaves as an ideal gas of 
clans, an interpretation which fits very nicely with the 
idea that clans arc independent objects, as implied by 
the definition of CPDs. 

The NB(Pascal)MD, with parameters n and k: 



p{n) 



fc(fc + l)...(fc + n-l) 



(14) 



is an example of a PSD, with the following identification: 
_k{k + l)...{k + n~l) 



(15) 



7(6) = {l-b) 



-k 
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Furthermore, the NBMD also belongs to the class of dis- 
crete IDD, the multiplicity distribution inside each clan 
being of logarithmic type. 

We obtain therefore the following value for N: 



N = -fcln(l - 6); 



(16) 



which also gives the the grand-canonical partition func- 
tion, applying eq. (pi): 



Q-(l-6)- 



(17) 



Comparing now with our proposed correspondence, 
eq. (p^), we find that fc is oi, i.e., the canonical partition 
function for a system with 1 particle, and must therefore 
be function of V and T: k — k{V, T); we also find that b 
is the fugacity of the system, i.e., b = exp(/i/fcBT), and 
it is a scaling function of n/k, see eq. (p^. 

We calculate now, using the standard thermodynami- 
cal relations, the average number of particles in the sys- 
tem, (n) , which turns out to be equal to the n parameter 
of the NBMD: 



keT 
ksT 



k b 
1 -bhiTr 1 - b 



-ksTk 
kb 



91n(l - b) 



T,V 



(18) 



consistently with the above mentioned relation b — 
n/{n + k). We also obtain that the average number of 
particles per clan, fic, is a function only of the fugacity 
of the system: 



(19) 



(fe- l)ln(l -6) 



This result is very interesting because formally, as already 
pointed out, 6 is a scaling function of n/k, and experi- 
mentally fic is seen to vary with the width of the rapid- 
ity interval at fixed cm. energy: if one were to identify 
(pseudo-)rapidity with volume, one would expect that 
changing At; at fixed cm. energy would imply changing 
the volume at constant temperature and fugacity (inten- 
sive variables), thus keeping fic constant, contrary to ob- 
servations. We must conclude that in the present ap- 
proach we cannot identify rapidity with volume as a sim- 
ple thought would suggest but we should allow volume 
to vary also with other physical quantities. 

We now turn our attention to the generating function 
(g.f.) for the multiplicity distribution, defined in eq. (|^). 
In the general case illustrated by eq. (|^) we easily find 



G{u) 



Q{uz,V,T) _ p(0)| 



Q(z,y,T) p(0)|„. 



(20) 



which is the ratio of the grand-canonical partition func- 
tion (or of the void probability) at two different fugaci- 
ties. This is a very general expression valid for any system 
in the grand-canonical ensemble. 



The g.f. for a CPD can always be written as 

GcPD(u) = exp[7Vg(?/)-iV] , (21) 

where g{u) is the g.f. for the multiplicity distribution 
within each clan (it satisfies g{Q) — 0). 

Because of eq. (E2|) we can write for the class of CPDs: 



G(u) =e-^^Q(uz,X^,r). 



(22) 



However, In Q{uz, V, T) — N{uz, V, T), hence we obtain 

G{u) = exp [N{uz) - N{z)] (23) 

which can be interpreted as a function of the difference in 
the average numbers of clans for a system with fugacity 
uz to that for a system at the actual fugacity z, keeping 
the same volume and temperature. For the g.f. within 
one clan we further find: 



N{uz) _ n{uz,V,T) _ P{uz,V,T) 
N{z) ~ n{z,V,T) ~ P{z,V,T) ■ 



(24) 



Interestingly, this is the ratio of the average number of 
clans for two systems with unequal fugacities. 

In addition it is interesting to remark that, remem- 
bering that parameter k depends on V and T, a com- 
plete thermodynamics can be built in the just mentioned 
framework. Its main quantities are listed in the following 
and explicitly calculated in the Appendix. 

The equation of state is 



PV 
k^ 



fcln 1 



(n) 
k 



(25) 



thus the average number of clans can be expressed in 
terms of the thermodynamic potential f2: 



N = -n/ksT. 



(26) 



The Helmholtz free energy can be rewritten in a form 
symmetric in h and k: 



A 



{n) In 1 + 



ksT \ (n) 

The average internal energy is 

U ,-/51nfc 



fcln 1 



knT 



N 



d\nT 



(27) 



(28) 



The entropy is 



S = kH<- 



A 
k^ 



(n) 
k 



V • (29) 



which coincides with —A/T in the limit of {dk/dT)v — > 
0, which gives also C/ ^ 0. 

In the next section we will focus our attention on clan 
thermodynamics of final charged particle MD. 
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IV. CLAN THERMODYNAMICS AND THE 
NB(PASCAL)MD 

In this Section an attempt is made in order to interpret 
in the present approach a surprising finding in some of the 
possible scenarios for hadron hadron coUisions in the TeV 
region discussed in references [Q, ^, i.e. the unexpected 
decrease with cm. energy of the average number of clans 
for semihard events in scenarios 2 and 3. 

We are guided by two considerations. 1) the occurrence 
of the NB(Pascal)MD — as it is the case in the scenarios 
mentioned above both for semihard and soft events — is 
usually interpreted as the result of a two step process: 
to the independent production of clans during the first 
step, it follows their decay according to a logarithmic dis- 
tribution, which can be obtained by a weighted average 
of geometric (Bose-Einstein) distributions during the sec- 
ond step. 2) the validity of the generalised local parton 
hadron duality (GLPHD). 

It should be pointed out that clan ancestors are in- 
dependently produced and Poissonianly distributed, by 
assumption, and a clan is, by definition, a group of par- 
tons of common ancestor; a clan consists of at least one 
parton, its ancestor. Each ancestor can be considered 
as an independent intermediate gluon source. All cor- 
relations among generated partons are exhausted within 
each clan. 

Clan ancestors can be produced either very early in 
the production process at higher virtualities or later at 
lower virtualities . 

In the first case, the ancestor's "temperature" (an 
unknown function in this approach of the average px 
and the rapidity) is expected to be higher: this ex- 
pectation, together with the lack of mutual correlations 
among ancestors, emphasises their overall quasi-classical 
behaviour: ancestor production in this case is compet- 
itive with the increase of gluon population within each 
clan. This situation is qualitatively closer to that ex- 
pected at hadron level for soft events and semihard events 
in scenario 1. 

In the second case ancestors are produced later, at 
lower virtualities: their "temperature" should also be 
lower, with even "colder" generated gluons. Their vir- 
tuality is lower. Accordingly, quantum effects should be 
expected to be enhanced in events sharing these proper- 
ties: new produced gluons prefer to stay together with 
other relatives within each clan than to become ancestor 
and initiate a new clan, k parameter is in this case lower 
and closer to that of a Bose-Einstein distribution, which 
occurs for k = 1. These remarks are consistent with the 
interpretation of (see Q ) as a measure of aggregation 
of partons into clans: it corresponds to the ratio of the 
probability to have two gluons (particles at hadron level) 
in the same clan over the probability to have two gluons 
in two separate clans, i.e., to smaller k parameter cor- 
responds an higher aggregation among produced gluons 
into clans. In addition, being 1/k linked to the integral 
of two parton rapidity correlations via second order fac- 



torial cumulants, the decrease of k implies stronger two 
parton correlations. In conclusion generated gluons pre- 
fer to stay together than to stay far apart, higher parton 
density regions are generated, the probability to create a 
new gluon is enhanced (a typical quantum effect), clans 
become more populated and their average number is re- 
duced. This situation is closer qualitatively to that ex- 
pected at hadron level for semihard events in scenarios 2 
and 3 of references [0, |j . 

The just mentioned considerations and Eq.s (p^, ( |T7| ) 
and (|l^) of Section |l| fully outline in the present ap- 
proach the importance of the behaviour of the fugacity 
variable b. 

The remaining question is how clan thermodynamics 
results at parton level can be extended to final particles 
through the hadronization mechanism. A possible an- 
swer to this question comes from generalised local parton- 
hadron duality (GLPHD) which says that all inclu- 
sive distributions are proportional at the two levels of 
investigation: 

Qn.hadrons (yi 1 ■ ■ ■ i l/ri) — P Qn. partons {yi,---,yn), (30) 
which corresponds for NBMD parameters to 

^hadron — ^parton: ^hadron — P^parton- (31) 

GLPHD can be applied separately to soft and semihard 
components thus solving our problem. In particular in 
this framework minijets production is related to the exis- 
tence of regions of high gluon densities and final particle 
production should be sensitive to the mentioned quan- 
tum effects, by increasing two particle correlations and 
BE effects. 

Motivated by these considerations, the behaviour of 
parameter & as a function of cm. energy, as well as the b 
dependence of N and fic have been explored as a sugges- 
tive example in the different above mentioned scenarios 
phenomenologically described in terms of NBMD's. In 
addition in view of their simple connections with N and 
fic, the probability of having no particles in the event, 
p(0), and the void scaling function, V{n/k), have been 
studied as a function of fugacity b. It turns out that the 
analysis in terms of p{0) and V{n/k) variables confirms 
the main result of the new approach, i.e., that the reduc- 
tion of the average number of clans with the increase of 
cm. energy is a quantum effect. 

We proceed now to discuss the thermodynamical be- 
haviour of multiparticle production according to the sce- 
narios described in the introduction and fully charac- 
terised by the cm. energy dependence shown in Fig. ^. 
We start with the soft component thermodynamical be- 
haviour since it is assumed to be the same in all above 
mentioned scenarios (see Fig. 

We observe that fugacity b is growing very fast from 
0.25 to 0.75 with cm. energies below 100 GeV and then 
smoothly varying from 0.75 up to 0.9 at 14 TeV. The 
explanation of this behaviour will be given in terms of 
the following three known possible interpretation of the 
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FIG. 1: Cm. energy dependence of standard NBMD parameters n and k in the three scenarios described in the introduction. 
The top two panes show the behaviour of n (it is the same in all scenarios for both the semihard and the soft component). 
The lower four panes show the k parameter for the soft component (equal in all scenarios) and the semihard one in the three 
scenarios. 
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FIG. 2; Results for the soft component (equal in all scenarios). In this figure, the lines and the open points show the results 
from our extrapolations: the solid line refers to full phase space, the dashed lines and open square to the interval |7;| < 1. The 
last point on each line correspond to a cm. energy of 14 TeV. The solid circles show full phase space data from ISR and UA5, 
the solid triangles refer to UA5 data in the interval \ri\ < 1. (a) fugacity fe as a function of cm. energy; (b) a parameter as a 
function of cm. energy; (c) average number of clans vs fugacity; notice that this is also a plot of the grand partition function 
in logarithmic scale, since N = log Q, Eq. (^^; (d) average number of particles per clans vs fugacity; here for clarity the curves 
and the data for \ri\ < 1 are shifted up by 2 units; (e) void probability vs fugacity; (f) void scaling function V vs fugacity; also 
here curves and the data for I77I < 1 are shifted down by 0.2 units. 
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FIG. 3: Same as figure 0, but the semihard component in scenario 1 is shown. 
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FIG. 4: Same as figure 0, but the semihard component in scenario 2 is shown. 



9 






0.8 1 

fugacity b 



FIG. 5: Same as figure 0, but the semihard component in scenario 3 is shown. 
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occurrence on NB (Pascal) MD in high energy physics phe- 
nomenology, which in terms of parameters n and k are: 



a) 



{n + l)p{n + l) 
p{n) 



a + [3n 



(32a) 



where a = kn/{n + k) and /? = n/{n + k). Notice that 
for a = (3, i.e. for fc = 1, the MD p{n) becomes a Bose- 
Einstein distribution, for P — 0, i.e. k — > oo, a Poisso- 
nian distribution (the Poissonian limit) and for a = a 
logarithmic distribution, which can be expressed as the 
superposition of Bose-Einstein distributions. 



^-= fcln(l + n/fc) 

(32b) 

where N is the average number of clans and fic the aver- 
age number of particles per clan. 



p(0) = 



(32c) 



where p{Q) is the probability of generating zero charged 
particles and V(n/fc) is the void scaling function; the oc- 
currence of scaling in the product of the first two mo- 
ments n and 1/fc indicates two particle correlation dom- 
inance for hierarchical systems, and the distance from 
point one (Poisson limit V(0) = 1) on the scaling function 
is larger, more numerous and correlated are the part icles . 

It should be pointed out that (3 parameter in Eq. (32a) 
coincides with the fugacity b discussed above in our ther- 
modynamical approach and therefore 6, like V(n/fc), is a 
scaling function of n/k, and a parameter corresponds to 
the average charged multiplicity, n, for a classical system 
{k — > oo). 

In this sense the relative behaviour oi (3 — b an d a = kb 
as the cm. energy increases in view of Eq. (32a) and the 
discussion at the beginning of this Section, can be con- 
sidered an indication of the relative importance of a be- 
haviour closer to a quantum one, i.e. harder, with respect 
to a behaviour closer to a quasi-classical, i.e. softer, for 
a class of events. A very slow increase of b with cm. en- 
ergy and an almost constant behaviour oi a = kb is the 
main characteristic of the class of soft events as shown in 
Figures |^a and |^. 

This fact is confirmed by inspection of Fig. ^ and d, 
where it is shown that is a very slow growing function 
of the fugacity of the system throughout the ISR region 
and below that region (« 7), and then a quickly growing 
function of the same variable in the GeV region up to 
14 in the TeV region (14 TeV); fie as a function of the 
fugacity has a similar behaviour from « 1.5 to « 3. 

Accordingly, the probability of creating zero charged 
particles, p(0), is decreasing throughout the same regions 
from 10~^ to 10~^ (for b — 0.9 at 14 TeV cm. energy); in 



addition (Fig. |^) the void scaling function V{n/k) turns 
out to populate for larger values of fugacity variable sec- 
tions of the curve far from the Poissonian limit (V(0) = 1) 
showing a clear increase of two particle correlations in 
this region as expected for a hierarchical system. 

It should be noticed that in the soft scenario (no mini- 
jets) constant k parameter behaviour in rapidity intervals 
as requested by KNO scaling in the GeV and TeV regions 
implies that also the other variables remain constant in 
the same regions. 

In scenario 1, as already pointed out, the semihard 
component is assumed to have a very similar behaviour to 
the soft one: KNO scaling is satisfied and minor changes 
in the general trend of the variables both in full phase 
space and in (pseudo)rapidity intervals are straightfor- 
ward consequences of the smaller constant k parameter 
value suggested by NB fits for the semihard component 
in the GeV and extrapolated to the TeV region (Fig. ||) . 

Coming to the second scenario the assumption of 
strong KNO scaling violation for the semihard compo- 
nent (an extreme point of view with respect to that of 
scenario 1) implies a completely new panorama with dra- 
matic changes. Fugacity b (Fig. ^a) is growing very fast 
from 0.4 at 200 GeV cm. energy up to 0.96 at 14 TeV 
almost saturating the maximum allowed value, which is 
one, and a parameter (Fig. ^) is decreasing very rapidly 
from « 16 at 200 GeV to « 3 at 14 TeV. 

The combined information contained in the two figures 
leads to the same conclusion, i.e. the proposed semihard 
scenario behaviour is much closer to a quantum one than 
the soft scenario favouring the production of regions of 
higher particle density. This interpretation is confirmed 
by studying general trends of the other variables as a 
function of fugacity b. The average number of clans N 
is decreasing in full phase space from « 30 (6 « 0.35) to 
w 10 at 14 TeV (6 w 0.96) and the average number per 
clan, fic is increasing from w 1 to « 8 in the same interval. 
Accordingly the probability of zero particle production is 
increasing from « 10"^'^ to « 10~^, i.e. gap probability 
is increasing with fugacity and cm. energy; in parallel 
void scaling function V{n/k) is populating sections with 
much higher b values than in scenario 1 corresponding 
to regions much farther from the Poissonian limit. One 
interesting point concerns smaller (pseudo) rapidity inter- 
vals (say 1 77 1 < 1): the general trend is that N is lower 
than in full phase space as are corresponding fic values, 
thus suggesting the onset of regions with higher parti- 
cle densities and lower temperatures. The probability of 
generating zero particles, in view of the higher densities, 
is therefore much higher and the void probability is far 
from the Poisson limit. 

Scenario 3 is a QCD inspired scenario for the semihard 
component: it assumes for parameter k a QCD behaviour 
(see Fig. |^). This scenario gives a panorama for our vari- 
ables which is intermediate between the two extremes, 1 
and 2. Fugacity b is increasing very fast with cm. energy 
as in scenario 2, but a parameter has a sweeter trend (it 
is « 6 at 14 TeV) indicating stronger independent pro- 
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duction; this fact is clearly shown in Fig ||c,d where N 
is larger and smaller than in scenario 2. Differences 
in p(0) and V{n/k) behaviours in full phase space as well 
as in (pseudo) rapidity intervals with respect to scenario 
2 are all consequences of the just mentioned remarks. 



CONCLUSIONS 



APPENDIX A 

The main quantities of clan thermodynamics are ex- 
plicitly calculated in the following. 
The Helmholtz free energy is 

A = {n)n -PV = j-^M - keTN. (Al) 



Clan thermodynamics has been investigated in order 
to explain the decrease with cm. energy of N for the 
semihard component and its increase for the soft one in 
the most realistic scenarios of multiparticle production in 
the GeV and TeV regions in hh collisions. It turns out 
that these two behaviours for clans point out structures 
closer to classical or quantum properties of gluons. A 
thermodynamical approach to multiparticle production 
was constructed on this basis. Results were determined 
in the framework of NB(Pascal)MD applied separately in 
the two components and can be extended to any infinitely 
divisible distribution. 



The average internal energy is 



dT 



b,V 

2 1 dk\ 



(A2) 



The entropy 



is 



S^^^ksN-{n)l^ + k,T^(^) = 
^..« + -=.Wl„(l + i)+.eT(|i)^,n(l + M 



The specific heat at constant volume is 



1 



0^ T^(^^^\ 



V 



(n) 



k{k + (n)) ■ 



(A3) 



(A4) 
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